SPARTA (7 Jul 2014)
# advect particles via VSS collisional flow on a uniform grid
# particles reflect off global box boundaries

variable            x index 10
variable            y index 10
variable            z index 10

variable            lx equal $x*1.0e-5
variable            lx equal 40*1.0e-5
variable            ly equal $y*1.0e-5
variable            ly equal 50*1.0e-5
variable            lz equal $z*1.0e-5
variable            lz equal 50*1.0e-5

variable            n equal 10*$x*$y*$z
variable            n equal 10*40*$y*$z
variable            n equal 10*40*50*$z
variable            n equal 10*40*50*50

seed	    	    12345
dimension   	    3
global              gridcut 1.0e-5

boundary	    rr rr rr

create_box  	    0 ${lx} 0 ${ly} 0 ${lz}
create_box  	    0 0.0004 0 ${ly} 0 ${lz}
create_box  	    0 0.0004 0 0.0005 0 ${lz}
create_box  	    0 0.0004 0 0.0005 0 0.0005
Created orthogonal box = (0 0 0) to (0.0004 0.0005 0.0005)
create_grid 	    $x $y $z
create_grid 	    40 $y $z
create_grid 	    40 50 $z
create_grid 	    40 50 50
Created 100000 child grid cells
  parent cells = 1
  CPU time = 0.0161159 secs
  create/ghost percent = 83.3893 16.6107

balance_grid        rcb part
Balance grid migrated 0 cells
  CPU time = 0.219159 secs
  reassign/sort/migrate/ghost percent = 4.74772 0.37967 0.751507 94.1211

species		    ar.species Ar
mixture		    air Ar vstream 0.0 0.0 0.0 temp 273.15

global              nrho 7.07043E22
global              fnum 7.07043E6

collide		    vss air ar.vss

create_particles    air n $n
create_particles    air n 1000000
Created 1000000 particles
  CPU time = 0.192408 secs

stats		    10
compute             temp temp
stats_style	    step cpu np nattempt ncoll c_temp

# first equilibrate with large timestep to unsort particles
# then benchmark with normal timestep

timestep 	    7.00E-8
run                 30
Memory usage per proc in Mbytes:
  particles (ave,min,max) = 102.997 102.997 102.997
  grid      (ave,min,max) = 17.624 17.624 17.624
  surf      (ave,min,max) = 0 0 0
  total     (ave,min,max) = 120.621 120.621 120.621
Step CPU Np Natt Ncoll temp 
       0            0  1000000        0        0    273.28433 
      10    4.6545749  1000000   963874   704709    273.28433 
      20    9.9629788  1000000   990032   705925    273.28433 
      30    15.551347  1000000  1003189   706182    273.28433 
Loop time of 15.5514 on 1 procs for 30 steps with 1000000 particles

Particle moves    = 30000000 (30M)
Cells touched     = 147311369 (147M)
Particle comms    = 0 (0K)
Boundary collides = 2597328 (2.6M)
Boundary exits    = 0 (0K)
SurfColl checks   = 0 (0K)
SurfColl occurs   = 0 (0K)
Collide attempts  = 29048355 (29M)
Collide occurs    = 21091819 (21.1M)
Particles stuck   = 0

Particle-moves/CPUsec/proc: 1.92909e+06
Particle-moves/step: 1e+06
Cell-touches/particle/step: 4.91038
Particle comm iterations/step: 1
Particle fraction communicated: 0
Particle fraction colliding with boundary: 0.0865776
Particle fraction exiting boundary: 0
Surface-checks/particle/step: 0
Surface-collisions/particle/step: 0
Collision-attempts/particle/step: 0.968279
Collisions/particle/step: 0.703061

Move  time (%) = 7.68697 (49.4295)
Coll  time (%) = 7.05109 (45.3406)
Sort  time (%) = 0.782014 (5.02858)
Comm  time (%) = 0.000120401 (0.000774217)
Outpt time (%) = 0.0311291 (0.20017)
Other time (%) = 5.50747e-05 (0.000354147)

Particles: 1e+06 ave 1e+06 max 1e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Cells:      100000 ave 100000 max 100000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
GhostCell: 0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
EmptyCell: 0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
timestep 	    7.00E-9
run 		    100
Memory usage per proc in Mbytes:
  particles (ave,min,max) = 102.997 102.997 102.997
  grid      (ave,min,max) = 17.624 17.624 17.624
  surf      (ave,min,max) = 0 0 0
  total     (ave,min,max) = 120.621 120.621 120.621
Step CPU Np Natt Ncoll temp 
      30            0  1000000  1003189   706182    273.28433 
      40    1.4746571  1000000    90285    70462    273.28433 
      50    2.9732702  1000000    91935    70649    273.28433 
      60    4.4813762  1000000    92973    70574    273.28433 
      70    5.9920361  1000000    94099    70726    273.28433 
      80     7.510874  1000000    94776    70401    273.28433 
      90     9.035507  1000000    95373    70527    273.28433 
     100    10.566248  1000000    96008    70701    273.28433 
     110    12.097086  1000000    96363    70715    273.28433 
     120    13.634879  1000000    96797    70549    273.28433 
     130    15.178449  1000000    97568    70886    273.28433 
Loop time of 15.1785 on 1 procs for 100 steps with 1000000 particles

Particle moves    = 100000000 (100M)
Cells touched     = 139109366 (139M)
Particle comms    = 0 (0K)
Boundary collides = 865894 (0.866M)
Boundary exits    = 0 (0K)
SurfColl checks   = 0 (0K)
SurfColl occurs   = 0 (0K)
Collide attempts  = 9371431 (9.37M)
Collide occurs    = 7021537 (7.02M)
Particles stuck   = 0

Particle-moves/CPUsec/proc: 6.58828e+06
Particle-moves/step: 1e+06
Cell-touches/particle/step: 1.39109
Particle comm iterations/step: 1
Particle fraction communicated: 0
Particle fraction colliding with boundary: 0.00865894
Particle fraction exiting boundary: 0
Surface-checks/particle/step: 0
Surface-collisions/particle/step: 0
Collision-attempts/particle/step: 0.0937143
Collisions/particle/step: 0.0702154

Move  time (%) = 7.43831 (49.0057)
Coll  time (%) = 4.93653 (32.5232)
Sort  time (%) = 2.69946 (17.7848)
Comm  time (%) = 0.000358105 (0.00235929)
Outpt time (%) = 0.103626 (0.682717)
Other time (%) = 0.000186682 (0.00122991)

Particles: 1e+06 ave 1e+06 max 1e+06 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Cells:      100000 ave 100000 max 100000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
GhostCell: 0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
EmptyCell: 0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0

